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Abstract. Aerosol perturbations over selected land regions 
are imposed in Version-4 of the Goddard Earth Observing 
System (GEOS-4) general circulation model (GCM) to as- 
sess the influence of increasing aerosol concentrations on re- 
gional circulation patterns and precipitation in four selected 
regions: India, Africa, and North and South America. Part 1 
of this paper addresses the responses to aerosol perturbations 
in India and Africa. This paper presents the same for aerosol 
perturbations over the Americas. GEOS-4 is forced with pre- 
scribed aerosols based on climatological data, which inter- 
act with clouds using a prognostic scheme for cloud micro- 
physics including aerosol nucleation of water and ice cloud 
hydrometeors. In clear-sky conditions the aerosols inter- 
act with radiation. Thus the model includes comprehensive 
physics describing the aerosol direct and indirect effects on 
climate (hereafter ADE and AIE respectively). Each simula- 
tion is started from analyzed initial conditions for 1 May and 
was integrated through June-July-August of each of the six 
years: 1982-1987 to provide a 6-ensemble set. Results are 
presented for the difference between simulations with dou- 
ble the climatological aerosol concentration and one-half the 
climatological aerosol concentration for three experiments: 
two where the ADE and AIE are applied separately and one 
in which both the ADE and AIE are applied. The ADE and 
AIE both yield reductions in net radiation at the top of the 
atmosphere and surface while the direct absorption of short- 
wave radiation contributes a net radiative heating in the atmo- 
sphere. A large net heating of the atmosphere is also apparent 
over the subtropical North Atlantic Ocean that is attributable 
to the large aerosol perturbation imposed over Africa. This 
atmospheric warming and the depression of the surface pres- 
sure over North America contribute to a northward shift of 
the inter-Tropical Convergence Zone over northern South 
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America, an increase in precipitation over Central America 
and the Caribbean, and an enhancement of convergence in 
the North American monsoon region. 

Keywords. Atmospheric composition and structure 
(Aerosols and particles; Cloud physics and chemistry) - Me- 
teorology and atmospheric dynamics (General circulation) 


1 Introduction 

The radiative impacts of aerosols can affect regional circula- 
tion patterns and their related hydrological cycling. Chung 
et al. (2002) impose the observed winter season aerosol ra- 
diative forcing on the Indian Ocean in a global model and 
report a reduction in winter precipitation over India. During 
the spring and early summer the clear-sky radiative forcing of 
dust and soot aerosols provides an elevated heat source over 
the Indo-Gangetic Plain of India that promotes the onset of 
the summer monsoon and increases precipitation in North- 
ern India (Lau and Kim, 2006; Lau et al., 2006). Yoshioka et 
al. (2007) attribute a portion of reduced precipitation along 
the inter-Tropical Convergence Zone (ITCZ) over Africa to 
the clear-sky aerosol radiative forcing, while Rotstayn and 
Lohman (2002) report a reduction in Sahel precipitation in 
response to the cloud radiative effect resulting from the en- 
trainment of pollution over the North Atlantic Ocean and the 
accompanying modification of cloud cover and cloud micro- 
physical processes. 

Little is yet known regarding the response over North and 
South America to the radiative impacts of aerosols. The 
smoke that accompanies forest and agricultural fires during 
the austral winter dry season is significantly absorbing (sin- 
gle scatter albedo approximately 0.9, Dubovik et al., 2002) 
and act as efficient cloud condensation nuclei (CCN; Hobbs 
andRadke, 1969; Kaufman and Fraser, 1997). The absorbing 
nature of biomass burning aerosols contributes a significant 
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Table 1 . Synopsis of simulation experiments. 


Experiment# 

Aerosol anomaly in selected regions: 
2X or IX or 1/2X of 
climatological values 

Aerosol Direct 
Effect (ADE) 
On/Off 

Aerosol Indirect 
Effect (AIE) 
On/Off 

1 

1 X 

on 

on 

2a 

2 X 

on 

on 

2b 

V 2 X 

on 

on 

3a 

2 X 

on 

off 

3b 

V2X 

on 

off 

4a 

2 X 

off 

on 

4b 

V 2 X 

off 

on 


Regional Mask 




PBL CCN {§ mm -3 ) MIRAGE 



1 SOW 1 20W 90W 60W 30W 0 30E 60E 90E 120E 150E 


Fig. 1 . Regions over which aerosol amounts are perturbed (top), 
Column aerosol optical depth from GOCART climatology (middle; 
Chin et al., 2002) used in ADE calculation, boundary layer CCN 
concentration from the MIRAGE climatology (bottom; Easter et al., 
2004) used in the AIE calculation. 


atmospheric heating and surface cooling which should in- 
duce a significant response in the surface energy balance and 
lower tropospheric stability. This can potentially yield sig- 


nificant increases or decreases in cloud cover depending on 
the level of smoke loading (Koren et al., 2008). Zhang et 
al. (2008) use a regional climate model to simulate the com- 
bined direct and indirect effects of smoke over South Amer- 
ica and find that smoke, and dynamical changes in cloud 
cover attributable to the smoke, delays the onset of the wet 
season and enhances moisture convergence and cloud cover 
upwind of the smoke region at the expense of the local mois- 
ture convergence. 

Under clear-sky conditions, aerosols scatter and absorb so- 
lar radiation. The scattering effect cools the column by ap- 
proximately similar magnitudes at the top of the atmosphere 
(TOA) and surface (Haywood et al., 1999). The absorbing 
effect yields a warming in the atmospheric layer at which 
the absorbing aerosol resides and yields a cooling of nearly 
equal magnitude at the surface (Haywood et al., 1999; Ra- 
manathan et al., 2001). These scattering and absorption ef- 
fects of aerosols are typically referred to as the aerosol direct 
effect (ADE). Aerosols entrained into clouds can modify the 
microphysical properties of the cloud if the aerosol particles 
act as CCN. A host of radiative consequences have been pro- 
posed to result from this interaction, commonly referred to 
as the indirect effects or (AIE). Assuming the same amount 
of condensed water is created by cloud dynamics, an envi- 
ronment that activates a greater number of CCN will have 
smaller drops and a greater albedo (Twomey, 1974; Hobbs, 
1999; Wilcox et al., 2006; and others). Smaller drop sizes 
inhibits the collision coalescence of liquid drops and thereby 
inhibits rainfall. This can lead to a radiative impact by in- 
creasing the lifetime of the cloud (Albrecht, 1989). An ad- 
ditional process by which aerosols motivate a cloud radiative 
impact is the semi-direct effect, whereby the ADE changes 
the stability of a cloud layer, resulting in a change (usually a 
reduction) in cloud cover (Hansen et al., 1997). 

The net effect of aerosols on climate may result from 
complex nonlinear interactions between the aerosol radia- 
tive perturbation, clouds, and the regional circulation. For 
example, where aerosol processes affect the vertical tem- 
perature profile of the atmosphere, that in turn influences 
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Table 2. Energy and water flux differences 2X minus 1/2X climatological aerosol amount for North and South America. 


Difference fields (2X minus 

ADE & AIE (exp. 2) 

ADE only (exp. 3) 

AIE only (exp. 4) 

Nonlinear interactions by difference 

1/2X) aerosol simulations 

N. Amer./S. Amer. 

N. Amer./S. Amer. 

N. Amer./S. Amer. 

N. Amer./S. Amer. 

Precipitation (mm/mon) 

q 

d 

—3. 7/1.0 

— 1.0/— 3.2 

4.8/— 1.8 

Evaporation (mm/mon) 

— 1.5/— 4.2 

—2.2/— 0.7 

—2.1/— 3.7 

2. 8/0.2 

Moisture Convergence 

(mm/mon) 

1.4/— 0.1 

— 1.5/1. 7 

0.8/0. 1 

2.1/— 1.9 

Sensible heat Flux 

—9.0/— 8.1 

—4.0/— 5.6 

—4.5/— 3.0 

— 0.5/0.5 

(W m -2 ) 

Latent heat Flux (W m~ 2 ) 

— 1.3/— 3.7 

— 1.9/— 0.6 

— 1.8/— 3.2 

2.4/0. 1 

TOA Shortwave (W m -2 ) 

— 10.5/— 10.1 

— 1.2/— 2.1 

— 10.7/— 9.3 

1.4/1. 3 

TOA net longwave (out- 
going minus incoming; 

— 1.8/— 2.5 

2.3/0.2 

—4.9/— 4.2 

0.8/1. 5 

W m -2 ) 

TOA Net Radtn. (Wm“ 2 ) 

—8.8/— 7.6 

—3.5/— 2.3 

—5.8/— 5.1 

0.5/— 0.2 

Surface Shortwave 

— 15.9/— 15.1 

—5.8/— 8.6 

— 11.0/— 9.3 

0.9/2. 8 

(W m -2 ) 

Surface net longwave 

—5.1/— 3.5 

—0.1/— 2.1 

—4.0/— 2.8 

— 1.0/1. 4 

(outgoing minus incoming; 
W m -2 ) 





Surface Net Radtn. 

— 10.8/— 12.6 

—5.7/— 6.6 

—7.0/— 6.5 

1. 9/0.5 

(W m -2 ) 






the thermodynamics of clouds and onset of moist convec- 
tion. These responses modulate vertical distributions of mass 
and moisture convergence, further altering the distribution of 
convection and cloud cover. An illustration of such responses 
is provided by Lau et al. (2006) who demonstrate how the 
ADE can strengthen monsoonal circulations through a feed- 
back in the atmospheric water cycle. 

While empirical evidence for these effects is mounting, 
complex nonlinear responses suggest that a complete under- 
standing of the role of these effects in global and regional 
climate requires a comprehensive numerical model with pa- 
rameterized physics for the nucleation of liquid and ice cloud 
particles, bulk cloud microphyics, and interactions between 
shortwave (SW) and longwave (LW) radiation in the atmo- 
sphere with aerosol and cloud particles. Towards this end, 
Sud and Lee (2007) have provided the implementation and 
evaluation of a model for aerosol nucleation of clouds for 
GCMs. Other efforts are ongoing at climate modeling cen- 
ters to evaluate the regional and global impacts of aerosols 
(e.g. Ming et al., 2007; Takemura et al., 2005; Lohmann et 
al., 2007). This paper represents a first set of sensitivity stud- 
ies designed to investigate the response of climate in sev- 
eral regions to perturbations in the aerosol loading of those 
regions. Simulations are performed where the summertime 
climatology of aerosols is halved and then doubled in four 
regions simultaneously: India, Africa, and North and South 
America. Part 1 of this paper (Sud et al., 2009) addresses the 
response to the aerosol perturbations in Africa and India. In 
Part 2 we address the same for North and South America. 


2 GEOS4 GCM and McRAS-AC aerosol-cloud param- 
eterization 

This study employs version-4 of the Goddard Earth Observ- 
ing System atmospheric general circulation model. The sim- 
ulations were performed using a horizontal resolution of 2.0 
degrees latitude x 2.5 longitude, and a total of 55 levels in 
the vertical with hybrid co-ordinates. A comprehensive AIE 
parameterization with consistent precipitation microphysics 
has been implemented in GEOS-4 and evaluated in single- 
column model simulations of the seasonal cycle of cloud 
cover at the Dept, of Energy Atmospheric Radiation Mea- 
surement Program instrumented site in the United States 
Southern Great Plains (Sud and Lee, 2007). The AIE sub- 
routines consist of four sub-models comprising i) the Foun- 
toukis and Nenes (2005) scheme for aerosol activation of liq- 
uid clouds, ii) the Liu and Penner (2005) scheme for het- 
erogeneous and homogeneous nucleation for ice-cloud parti- 
cles, iii) the Sud and Lee (2007) precipitation microphysics 
adapted from Seifert and Beheng (2001) for a coarse resolu- 
tion GCM, and iv) the Khvorostyanov and Curry (1999a, b) 
parameterization of the cloud particle effective radius. The 
AIE parameterization solves for cloud particle number den- 
sity and effective radius. These parameterizations are inte- 
grated into the existing Microphysics of Clouds with Re- 
laxed Arakawa- Schubert (McRAS) scheme for bulk cloud 
physics and microphysics in coarse-grid GCMs (Sud and 
Walker, 1999a, b), which determines cloud fraction and 
cloud condensate mass concentration based on the convective 
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Fig. 2. Difference in net shortwave radiation at surface (left) and top of the atmosphere (right) for doubled minus halved aerosol perturbations 
in Wm“ 2 . Top panels are for simulations with ADE only, middle panels are for simulations with AIE only and bottom panels are for 
simulations with ADE and AIE combined. 


and large-scale dynamics. The combined McRAS and AIE 
parameterization is referred to as McRAS -AC. The AIE is 
computed using the climatology of aerosol particle con- 
centrations from the MIRAGE model (Easter et al., 2004). 
McRAS- AC currently uses five aerosol species: black car- 
bon, sulfate, dust, organic carbon, and sea salt, but its algo- 
rithms are general enough to use any number of aerosols as 
long as their chemical properties are specified. 

The radiative transfer in GEOS-4 is described in Chou and 
Suarez (1994, 1999). The ADE is included using the optical 
properties of aerosols prescribed from the monthly aerosol 
data simulated with the Georgia Institute of Technology- 
Goddard Global Ozone Chemistry Aerosol Radiation and 
Transport (GOCART) model (Chin et al., 2002). These prop- 
erties assume that the aerosol species are externally mixed. 

Further details of the model and the relevant parameteriza- 
tions can be found in Part 1 of this paper (Sud et al., 2009), 
Sud and Lee (2007) and the references noted above. 


3 Simulation experiments 

This study focuses on the boreal summer season, which for 
South America is the dry winter season. Over North Amer- 
ica we have selected the region from the Great Plains to the 
East Coast, which is dominated by moist convective cloud 
systems and also experiences a rise in aerosol optical thick- 
ness in summer (Holben et al., 2001). The dry season over 
South America exhibits a dramatic rise in smoke concentra- 
tions as a consequence of extensive forest burning and low 
rainfall. Under clear-sky conditions this smoke yields a sub- 
stantial direct warming of the atmosphere and cooling of the 
surface. However, the net radiative effect of the aerosols de- 
pends strongly on the response of clouds through both the 
AIE and the semi-direct effect. Koren et al. (2008) find 
that cloud cover over the Amazon increases with increasing 
smoke aerosol by suppressing coalescence and rain, but also 
decreases with increasing aerosol because of semi-direct ef- 
fect. Which process wins out, they conclude, depends criti- 
cally on the amount of aerosol and the amount of cloud cover. 
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Sfc Net LW difference (W m -2 ) TOA Net LW difference (W rrr 2 ) 


ADE only 


ADE only 



In addition to the climate response in these specific regions, 
the combined effect of doubling the aerosol concentrations in 
these two regions and Africa during the same season yields 
a strong response in the tropics between North and South 
America. This response is also the subject of this study. 

Note that while the model accounts for the advection of 
altered cloud fields, the aerosol perturbations are limited to 
the land regions, therefore there is no accounting for the ad- 
vection of increased aerosol concentrations in these experi- 
ments. In this sense these experiments are sensitivity tests to 
a radiative perturbation in limited regions that is consistent 
with an increase in aerosols in that region. Of particular in- 
terest in this study is the atmospheric circulation response to 
the radiative perturbation arising from the forcing attributed 
to the ADE and AIE in a region. 

A summary of the model simulation experiments per- 
formed in this study is shown in Table 1. Experiment 1 
is a baseline simulation with climatological aerosols. Ex- 
periments 2 through 4 comprise 6-member ensembles each 
of simulations where the aerosol amounts are doubled and 
halved in the aerosol perturbation regions indicated in Fig. 1. 


Both the ADE and the AIE are active in the model in exper- 
iment 2, while the AIE and ADE are individually turned off 
in experiments 3 and 4 respectively. The climate response 
to these aerosol perturbations is evaluated from the differ- 
ence between the ensemble mean results of the doubled and 
halved aerosol perturbation simulations. 

All simulations used initial conditions derived from the 
NCEP reanalysis (Kalnay et al., 1996) interpolated to the 
GEOS-4 GCM grid for 1 May of each of the six years: 1982 
through 1987. Each simulation was integrated to 31 August 
to provide output to be examined for June- July- August (JJA). 
The first month (May) is ignored for initial adjustment to the 
imposed aerosol anomalies and intrinsic model biases. The 
initial conditions are chosen so that the model is started from 
a balanced state appropriate for the season. 1 May from each 
of six different years is used as a convenient means of gener- 
ating a six-member ensemble with different initial conditions 
for each ensemble member. This removes some short-term 
weather and interannual variability from the results. The fo- 
cus of the study is on the short-term response of the atmo- 
sphere to differences in aerosol amount over a single season, 
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Sfc Net Radiation difference (W m“ 2 ) TOA Net Radiation difference (W m" 2 ) 



Table 3. Cloud cover differences 2X minus 1/2X climatological aerosol amount for North and South America (—0 denotes a decrease of less 
than 0.5%). 


Difference fields (2X minus 1/2X) aerosol simulations. 

ADE & AIE (exp. 2) 
N. Amer./S. Amer. 

ADE only (exp. 3) 
N. Amer./S. Amer 

AIE only (exp. 4) 
N. Amer./S. Amer 

Low cloud (below 700 hPa level) cover (%) 

1. 3/2.7 

-0/0.1 

1. 8/2.4 

mid-level cloud (between 700 hPa and 400 hPa level) cover (%) 

0.5/1.9 

-0/0.2 

1. 1/2.8 

high cloud (above 400 hPa level) cover (%) 

-1.0/-0 

-1.0/- 1.0 

0.9/0. 2 

total cloud cover (%) 

0.2/3. 2 

— 2.0/— 0 

2.7/4. 7 


therefore all simulations used prescribed sea surface tem- 
peratures based on the Reynolds and Smith (1994) data set 
obtained from the Hadley Center. Climate responses over 
longer time scales associated with changes in sea surface 
temperature attributable to aerosol radiative forcing are not 
included in these results. 


4 Results 

4.1 Shortwave (SW) radiative fluxes at the top of the at- 
mosphere (TOA) and surface 

The aerosol perturbations are limited to the land areas in- 
dicated in Fig. 1. These are the regions where the short- 
wave flux difference attributable to aerosol scattering and in- 
creased cloud reflectivity are greatest. The magnitudes of 
the differences in radiative fluxes and water and energy cycle 
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components between doubled and halved aerosol amounts in 
these regions is summarized in Table 2. The magnitude of 
the change in net SW radiation at TOA attributable to the 
ADE is — 1.2 Wm -2 and — 2.1 Wm -2 in N. America and 
S. America, respectively, while the magnitude of the ATE is 
substantially greater; — 10.7 Wm -2 and — 9.3 Wm -2 . Maps 
of the change in TOA and surface net shortwave radiation 
are shown in Fig. 2. The change in net SW at the TOA in the 
experiment with ADE and AIE combined is similar to that 
of the AIE only experiment. However, the surface net SW 
reveals that the direct effect of aerosols remains an impor- 
tant process over the Americas in the simulations where the 
ADE and AIE are combined. In the ADE only simulation the 
change in net S W at the surface exceeds the change in net SW 
at the TOA by more than 5 W m -2 in both regions owing to 
atmospheric absorption. The difference in the change in net 
surface and TOA SW in the combined ADE and AIE exper- 
iment is comparable in magnitude. Since cloud absorption 
of SW radiation is small owing to the high single scattering 
albedo of clouds, the change in direct absorption of short- 
wave radiation is attributable to aerosols in the combined 
ADE and AIE experiment. 

While the change in radiative fluxes is largest over the 
aerosol perturbation regions, the circulation response and the 
advection of altered cloud fields yields significant changes in 
the radiative fluxes over the oceans. The circulation response 
(discussed further below) to the ADE over North America is 
a weakening of the westerly winds over the mid-latitude At- 
lantic Ocean (see Fig. 8b discussed further below) and a con- 
comitant reduction in evaporation and cloud cover. Hence 
there is a substantial increase in net SW radiation over the 
North Atlantic Ocean directly east of North America. This 
response is not present in the AIE experiment and is re- 
stricted to a smaller portion of the North Atlantic Ocean in 
the combined ADE and AIE experiment. 

4.2 Longwave (LW) and net radiative fluxes at the TOA 
and surface 

Longwave (LW) radiative responses in the ADE only exper- 
iment (Fig. 3, shown as outgoing minus incoming LW radi- 
ation) are largely attributable to changes in temperature and 
clouds in response to the aerosol SW radiative effect. Sulfate 
and smoke aerosol particles of North and South America do 
not significantly absorb in the longwave (LW) owing to their 
small size. At the TOA, the net LW change in the ADE only 
experiment is negligible in South America and 2.3 Wm -2 
in North America. This TOA cooling over North America 
is a consequence of the 2% decrease in cloud cover shown 
in Fig. 6 (percent changes in cloud cover are summarized 
in Table 3). In contrast, the change in LW at the surface 
in the ADE only experiment is negligible in North Amer- 
ica and — 2Wm -2 in South America. Unlike over North 
America, there is no net change in total cloud cover in the 
ADE only experiment over South America. The SW heat- 


Sfc Skin Temperature difference (C) 



Fig. 5. Same as Fig. 2 for surface skin temperature (in °C). No 
change is apparent over ocean because sea surface temperatures are 
fixed to observed values. 


ing owing to solar absorption by smoke in the atmosphere 
over South America exceeds 6 W m -2 in the ADE only ex- 
periment. This warming, together with a small increase in 
low cloud cover (<1%; see Fig. 7 and Table 3), enhances the 
downward LW radiation from the lower troposphere to the 
surface. Also, the decrease in surface temperature reduces 
the surface emission. Note, however, that the surface LW 
heating in South America is only one-fourth of the magni- 
tude of the SW cooling from smoke scattering and absorp- 
tion. The weak low cloud response over South America in 
the ADE only experiment suggests that the semi-direct effect 
observed in analyses of satellite data (e.g. Koren et al., 2008) 
is not driving a significant reduction in cloud cover in this 
model. 

LW heating is apparent at the TOA and surface in the 
AIE only simulations. In this case it is likely that the in- 
crease in CCN yields an increase in low and high cloud 
cover that both enhances the cloud greenhouse effect and 
enhances the downward LW radiation to the surface from 
cloud base. Net outgoing LW at the TOA is — 4.9Wm -2 
and —4.2 W m -2 over North America and South America re- 
spectively in the AIE only experiment. Surface net outgoing 
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High-level Cloud Fraction difference 


Column Cloud Fraction difference 


ADE only 


ADE only 



Fig. 6. Same as Fig. 2 for high cloud fraction (left) and column cloud fraction (right). High clouds are clouds at pressure levels less than 
400 hPa. 


LW is — 4.0Wm -2 and — 2.8 Wm -2 , respectively, indicat- 
ing that in response to the AIE LW heating over North Amer- 
ica is absorbed primarily by the surface, while over South 
America the net LW heating is divided roughly equally be- 
tween the surface and atmosphere. 

An increase in low cloud cover has been hypothesized in 
studies indicating a positive correlation of aerosol amounts 
and cloud cover assuming aerosol SW absorption is limited 
(Kaufman and Koren, 2006). Greater concentrations of CCN 
yield smaller drops for clouds of the same liquid water con- 
tent, and inhibit the collision coalescence of raindrops. This 
suppression of the warm rain process can increase the life- 
time, and hence mean coverage, of low clouds (Albrecht, 
1989). This process can also invigorate convection leading 
to more intense convection and an increase in the transport 
of liquid condensate to the detraining layers of convective 
clouds (Koren et al., 2005; Rosenfeld and Woodley, 2000). 
These processes, in addition to any increase in the trans- 
port of ice nucleating aerosols, are the primary microphys- 
ical causes of increases in cloud amount in simulations with 
increased aerosol amounts. However, it is important to note 


that the changes in cloud cover resulting in these simulations 
are also associated with changes in regional circulation pat- 
terns. Evidence of this is present in the substantial changes 
in LW and SW radiative fluxes over Central America where 
the aerosol concentrations were not changed, yet substantial 
changes in the circulation and precipitation are apparent (dis- 
cussed further below). 

In the simulations with combined ADE and AIE, net ra- 
diation at the TOA is a cooling over both North and South 
America (Fig. 4). The magnitudes are — 8.8 Wm -2 and 
— 7.6Wm -2 , respectively. Net radiation at the surface is a 
greater cooling of — 10.8 W m -2 and — 12.6 W m -2 for North 
and South America. The net heating of the atmosphere is 
smaller than that attributable to shortwave heating in the AIE 
only experiment reflecting both the modest net LW cooling of 
the atmosphere in both regions and the fact that the increase 
in cloud cover attributable to the AIE reduces the magnitude 
of the ADE when the two effects are combined. 

The net radiative cooling at the surface leads to a reduc- 
tion of the surface temperature (Fig. 5). Since the sea surface 
temperature is fixed to observed values in these simulations, 
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Mid— level Cloud Fraction difference 


Low-level Cloud Fraction difference 


ADE only 



ADE only 



0.2 


0.16 

0.12 



Fig. 7. Same as Fig. 2 for mid-level cloud fraction (left) and low cloud fraction (right). Low clouds are at pressure levels greater than 700 hPa. 
Mid-level clouds occur between 700 and 400 hPa. 


the reduction in surface temperature is restricted to the land 
areas. The cooling in the combined ADE and AIE experi- 
ment is — 0.67°C and — 0.86°C in North America and South 
America, respectively. 

In addition to the aerosol perturbations imposed over 
North and South America, aerosol perturbations are also im- 
posed over Africa and India in the same simulations. The 
perturbation in Saharan dust over North Africa yields an ex- 
tended response in the cloud radiative effect over the tropical 
North Atlantic Ocean attributable to the AIE in the AIE only 
simulations and the combined ADE and AIE simulations. 
This forcing is a likely contributor to the circulation and pre- 
cipitation response in Central America, the Caribbean, and 
northeastern South America. The net radiative effect of the 
AIE over the tropical North Atlantic Ocean is a heating of the 
atmosphere of greater than 15 Wm -2 . Moura and Shukla 
(1981) show that a heat source north of the equator and a 
heat sink south of the equator, in the form of sea surface tem- 
perature perturbations to the tropical Atlantic Ocean yields 
a northward displacement of the ITCZ that provides an in- 
crease in rainfall over the northernmost portion of South 


America and a reduction in precipitation in Northeastern 
Brazil. A similar pattern for the displacement of the ITCZ is 
apparent over the Atlantic Ocean and northern South Amer- 
ica in the AIE only and combined ADE and AIE simulations 
here (Fig. 9a discussed further below) which in this case is 
caused by the net atmospheric heating attributable to the LW 
component of the AIE over the Atlantic Ocean in response to 
a large aerosol perturbation over Africa. 

The origin of the radiative perturbation over the North At- 
lantic Ocean is the AIE over Africa. Although cloud cover 
North of the ITCZ over Africa is low, there is a 10% to 20% 
increase in high cloud cover over Africa resulting from the 
AIE. This occurs because dust acts as ice cloud nuclei in 
the McRAS-AC scheme. Strong mixing over the Sahara 
desert leads to a thick layer of dust in the aerosol clima- 
tology. Therefore, in our experiments we have halved and 
doubled the dust concentration leading to large differences 
in the number of ice nuclei in the troposphere. In the Sa- 
hel and continental ITCZ, additional ice nuclei may be con- 
tributed from the vertical transport of dust and other ice nu- 
clei in convective updrafts, which is also a process included 
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Fig. 8a. Difference in surface pressure for doubled minus halved 
aerosol perturbations. Top panel is for simulations with ADE only, 
middle panel is for simulations with AIE only and bottom panel is 
for simulations with ADE and AIE combined. 


Fig. 8b. Same as Fig. 8a for 10 m surface winds with shading indi 
eating convergence (in 10 -7 s -1 ). 


in the McRAS-AC scheme. Because of the additional ice nu- 
clei present when the dust concentration is increased, the ice 
water content of the cirrus clouds is composed of a greater 
number of smaller ice particles. The growth rate and settling 
rate of the ice particles both decrease with decreasing par- 
ticle size. This dependence on particle size means that the 
cirrus clouds in the simulations with higher dust concentra- 
tions persist longer, and mean cloud cover is greater. 

The McRAS cloud parameterization also allows for the 
advection of cloud condensate. The additional high cloud 
condensate attributable to the AIE over Africa is advected 
over the Atlantic Ocean owing to the upper- tropospheric 
winds. Thus, although there is no aerosol perturbation over 
the Atlantic Ocean in these simulations, the perturbation over 
North Africa induces a strong radiative perturbation over the 
ocean. The impact of cloud advection is generally small in 
GCMs relative to other components of cloud parameteriza- 
tions. Nevertheless, in a simulation of present-day climate 
the inclusion of this effect yields more realistic global cloud 
radiative effects (Sud and Walker, 1999b). Although it is dif- 


ficult at present to validate cloud advection in the model, we 
note that the erosion of cloud by mixing of dry air in the 
cloud layer, as well as evaporation driven by subsidence of 
the cloud layer or settling cloud particles, are processes in- 
cluded in McRAS-AC that act to limit the duration of ice 
cloud lifetime. We also note that a 4-times change in the 
magnitude of Saharan dust loading may be excessive com- 
pared to realistic scenarios. The halving and doubling of the 
aerosol amount in these experiments was done to perform 
a simple sensitivity study to aerosol perturbations of a sig- 
nificant and plausible magnitude compared to changes in in- 
dustrial emissions over the modern era. These perturbations 
were imposed on the dust in the interest of applying the ex- 
periment methodology consistently. To the extent that these 
dust perturbations may be excessive, so too is the response to 
the dust perturbation described here. 

4.3 Surface pressure and 10 m surface winds 

The change in surface pressure in the simulations with ADE 
and AIE combined is the depression of the North Ameri- 
can low and enhancement of the high-pressure cells over 
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Fig. 9a. Difference between doubled and halved aerosol amounts for precipitation (top), column moisture convergence (middle), and total 
precipitable water (bottom) for the ADE and AIE combined (left) and the AIE only (right). 
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Fig. 9b. Difference in precipitation between doubled and halved 
aerosol amounts for the ADE only. 


the Southeast Pacific and South Atlantic Oceans (Fig. 8a). 
These reinforce the seasonal patterns of surface pressure in- 
dicated in the long-term mean of the NCEP reanalysis (not 
shown). Over North America, the anomalous low pressure 
leads to a large southerly wind anomaly along 15 N (Fig. 8b). 
This wind anomaly amplifies the North American monsoon 
circulation in particular, and contributes to anomalous con- 


vergence over the American Southwest, as well as a large 
positive precipitation anomaly in Central America (Fig. 9a 
discussed further below). The subtropical meridional wind 
anomaly extends across the Americas from the East Pacific 
Ocean to the West Atlantic Ocean and is likely also respond- 
ing to the large thermal anomaly over the North Atlantic 
Ocean attributable to the AIE discussed above and implicated 
in the precipitation anomaly over South America. 

4.4 Precipitation, moisture convergence, and precip- 
itable water 

The change in precipitation, column moisture convergence, 
and precipitable water for doubled minus halved aerosols is 
shown for the ADE and AIE combined runs and the AIE only 
runs in Fig. 9a. The precipitation difference for the ADE 
only experiment is shown in Fig. 9b. The response of precip- 
itation to the ADE and AIE is a northward shift of the ITCZ 
and an enhancement of precipitation over Central America. 
The precipitation response to the perturbation in surface pres- 
sure over North America and the Atlantic Ocean is discussed 
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above. The circulation and associated precipitation perturba- 
tions are modest in the ADE only simulations and substan- 
tially greater in the AIE only and combined ADE and AIE 
simulations. Column moisture convergence anomalies co- 
incide with the precipitation anomalies reflecting the ITCZ 
shift and enhanced convection over Central America and the 
Caribbean. This leads to a large positive perturbation in pre- 
cipitable water throughout Central America, the Caribbean 
and the southern United States, and a decrease in precipitable 
water throughout western South America and the adjacent 
southeast Pacific Ocean. 


5 Summary 

The impact of aerosol perturbations on regional radiative 
forcing, circulation patterns, and precipitation is explored 
over the Americas during boreal summer in an ensemble 
set of GCM simulation experiments. The response is diag- 
nosed from the difference between simulations where cli- 
matological aerosol concentrations are doubled and halved 
over the selected land areas of North America, South Amer- 
ica, Africa, and India. The model includes parameterizations 
for radiative interactions of aerosols under clear- sky condi- 
tions (the ADE) and parameterized aerosol nucleation of liq- 
uid and ice clouds (the AIE). Separate experiments are per- 
formed to investigate the impacts of the ADE and AIE indi- 
vidually and combined. 

Both the ADE and the AIE yield net radiative cooling at 
the top of the atmosphere and the surface over both South 
America and North America. In both cases the surface forc- 
ing exceeds the TOA forcing, meaning there is a net radiative 
heating of the atmosphere. The heating results from direct 
aerosol absorption of solar radiation. The cloud greenhouse 
effect is enhanced by an increase in high cloud amount, how- 
ever the increase in low cloud amount attributable to the AIE 
leads to a greater LW heating of the surface and a net LW 
cooling of the atmosphere. At least a portion of the increase 
in high cloud amount, in particular the increase simulated 
over North Africa, is a consequence of an increase in ice nu- 
clei, as well as a suppression of the warm rain process in 
the lower part of convective clouds causing increased ver- 
tical transport of cloud condensate in convective updrafts. 
Changes in cloud amount also reflect changes in the regional 
circulation. In this study, we are unable to quantitatively sep- 
arate the contribution owing to the microphysical modifica- 
tion of clouds from any contribution owing to changes in the 
regional circulation. 

The response over the Americas is also likely impacted 
by the substantial response of clouds over the low latitude 
North Atlantic Ocean. A large LW heating perturbation re- 
sults from the advection of enhanced upper-level cloud gen- 
erated by the AIE over Africa in response to the increase in 
dust, which act as ice nuclei in the model. This cloud conden- 
sate is produced over the continent and transported over the 


Atlantic Ocean by the upper-tropo spheric winds. The LW 
heating perturbation over the ocean promotes a northward 
shift of the ITCZ over the Atlantic Ocean and northern South 
America. This feature, together with the depression of the 
seasonal low pressure minimum over North America results 
in a strong southerly wind anomaly that spans the Americas 
along 15 N latitude enhancing precipitation throughout Cen- 
tral America and enhancing convergence in the North Amer- 
ican monsoon region. 

These results represent a sensitivity study to strong aerosol 
perturbations in four selected regions of known high aerosol 
loading owing to industry, agriculture, and forest burning 
of anthropogenic origin. In light of the known long range 
transport of aerosols and the widespread response over sur- 
rounding oceans, and even the remote Atlantic Ocean, to 
aerosol perturbations over land, these experiments should be 
followed up with perturbations of global aerosols. Further- 
more, the results presented here may be dependent on the 
specific model physics implemented in the aerosol, cloud and 
convection parameterizations employed in this model. Addi- 
tional sensitivity tests of model parameters as well as inter- 
comparisons with other global models are warranted. 
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